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Abstract 

This work sets the non isotropic noncentral elliptical shape distributions via QR 
decomposition in the context of zonal polynomials, avoiding the invariant polyno- 
mials and the open problems for their computation. The new shape distributions 
are easily computable and then the inference procedure can be studied under ex- 
act densities instead under the published approximations and asymptotic densities 
under isotropic models. An application in Biology is studied under the classical 
gaussian approach and a two non gaussian models. 



1 Introduction and the main principle 

Considering the non isotropy in the non ce ntral shape theory ha s been very problematic, 
even in the gaussian case, see for example iGoodall and Mardia I (|1993l ) ; the correspond- 
ing shape densities involve expansions products of po wers o f trace s of different matrices, 
which forces the apparition of invariant polynomials ( Davis (jl908l )). Then the resulting 
densities enlarge the list of uncountable densities in the noncentral multivariate statistics, 
which can not be computable, and remains as theoretical results, very far from the infer- 
ence and the applications. So, the applications in shape theory have been force to avoid 
those polynomials, but at a very high cost, the assumption of isotropy. However, in this, 
the resulting densities we re series of zonal polynomia ls, and they could not studied prop- 
erly (before the works of Koev and Edelman ( 20061 )). and again t hey forced to the use 



of app r oximations and asymptotic distribution to perform inference, IGoodall and Mardia 
( 1993f ) . ISrvden and Mardia" ( 1998 ). and the references therein 



The following principle solves the first and more important problem, avoiding the in- 
variant polynomials, and setting the new shape distributions in terms of scries of zonal 
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polynomials, this series can be computed accurately and efficie ntly by simple modification s 
of the powerful algorithms of hypergeometric functions given bv lKoev and Edelman ( 2006 ). 



From the point of view of applications, the isotropic assumption = 1^ for an elliptical 
shape model of the form X ~ ^wxk(a* x5 S x , 0, h), restricts substantially the correlations 
of the landmarks in the figure. So, we expect the non isotropic model, with any positive 
definite matrix 0, as the best model for considering all the possible correlations among 
the anatomical (geometrical o mathematical) points. However, using the classical a pproach 
of the published literature of shape (see for example Goodall and Mardia ( 19931 )) under 



the non isotropic model, we obtain immediately invariant polynomials, which can not be 
computed at this time for large degrees. 

In order to avoid this problem consider the following procedure: Let 



X-5Ar xK (/X x ,S x ,0,/l), 

if 1 / 2 is the positive definite sq uare root of the matrix 0, i .e. = (0 1 / 2 ) 2 , with 1 



/2 • 



K x K, iGupta and Vargal (|1993l , p. 11), and noting that 



where 
then 



xe^x' = x(0- 1 / 2 ©- 1 / 2 )- 1 x' - x©- 1 / 2 (X0- 1 / 2 )' - zz', 

Z = X0~ 1/2 , 



Z ~ £ivxif (m z j S x ,Ix , h) 
with u = u !/ 2 , fsee lGupta and Vargal (|l993l p. 20)). 



And we arrive at the classical starting point in shape theory where the original landmark 
matrix is replaced by Z = X0~ 1//2 . Then we can proceed as usual, removing from Z, 
translation, scale, rotation and/or reflection in order to obtain the shape of Z (or X) via 
the QR decomposition, for example. 

Namely, the QR shape coordinates u of X are constructed in several steps summarized 
in the expression 

LX0 1/2 = LZ = Y = TH = rWH = rW(u)H, (1) 

which we discuss next. Observe that fi z = /x x */ 2 and the QR shape coordinates of fi z 
are defined analogously. The matrix L : (JV — 1) x N has orthonormal rows to 1 = (1, . . . , 1)'. 
L can be a submatrix of the Helmert matrix, for example. 

Let /x = L/x x , then Y : (N — 1) x K is invariant to translations of the figure Z, and 

Y ~ £ N - 1x k(v®~ 1/2 , S ® I K ,h), 

where £ = LS X L'. 

Now, let ben = min(iV— 1, K) andp = rank/i. In ([TJ) , Y = TH is the QR decomposition, 
where T : (JV— 1) x n is lower triangular with tu > 0, i = 1, . . . , min(n, K—l), and H : nxK, 
H € V n ,K , the Stiefel manifold. Note that T is invariant to translations and rotations of Z. 
The matrix T is referred as the QR size- and- shape and their elements are the QR size-and- 
shape coordinates of the original landmark data Z. Typically in shape analysis there are 
more landmarks than dimensions (N > K). H acts on the right to transform 3i K instead 
of acting on the left as in the multivariate analysis. In our case we see the landmarks as 
variables and the dimensions as observations, then the transposes of our matrices Z and Y 
can be seen as classical multivariate data matrices. 

According to the nature of the base H and providing that N — 1 > K, we say that 
T contains the QR reflection size- and- shape coordinates if H includes reflection, i.e. H g 
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0(K), |H| = ±1 and t KK > 0; otherwise, if H excludes reflection, H e SO{K), |H| = +1, 
tKK is not restricted, we say that T contains the QR size-and-shape coordinates. These 
cases will denote by T R and T iYfi , respectively. In the classical multivariate case, n < K, 
we do not have such classifications for T. 

Now, if we divide T by its size, the centroid size of Z, 



r = ||T|| = VtrT'T = j|Y||. 

we obtain the so called QR shape matrix W in (TTJ). Wc define W fl = T R /r or W"* = 
T NR /r if W includes or excludes reflection, respectively, and given that ||Wj| = 1, the 
elements of W are a direction vector for shape, and u comprises m = (N — 1)K — nK + 
\n(n + 1) — 1 generalized polar coordinates. 

Before deriving the main results of this paper we must solve some discrepancies between 
the shape theory and the classical multivariate theory. Recall that for a given Y : n x K, 
n = N — 1 < K, then YY' has the noncentral Wishart distribution which is invariant 
to orientation and reflection, but if n > K that density does not exist with respect to 
the Lebesgue measure defined on the space of positive definite n x n matrices, and we 
therefore use the size-and-shape matrix T. However, the density of YY' when, n > K, 
exist on the (nK — K(K — l)/2)-dimensional man ifold of rank-K positive semidefinit e n x n 



matri ces with K distinct positive eigenvalues, see iDiaz-Garcia and Gonzalez-Farias I (|2005f ) 



and Diaz-Garcfa and Gutierrez- Jaimez ( 20061 ). This last fact can provide an alternative 



form to study the shape theory, which is being analysed by the authors at present. 

And finally, classical integration over O(K) involving zonal polynomials gives the density 
of T R , but demands integration over SO(K), then we just recall that the corresponding 
integrals are the same when n < K, and, for n > K and p < K, they are twice the integral 
over SO{K). 

This work is distributed as follows: first, the size and shape distribution for any elliptical 
model with a full Kronecker covariance matrix is derived in section[2] The the shape density 
is obtained in section [H and the classical isotropic gaussian shape density, full derived in 



Goodall and Mardia I (|1993l ) , follows here as a corollary, then the section U describes the 



excluding reflection shape densities. The central case of the shape density is studied in 
section [5J and a remarkable property is established, i.e. it is established that the central 
QR reflection shape density is invariant under the elliptical family. Finally, some particular 
elliptical densities are derived in section [6] in order to perform inference on exact distribu- 
tions; a subfamily of Kotz distributions which contains the gaussian one is derived, then 
two elements of that class (the gaussian and a non gaussian model) is applied to an existing 
publish data, the mouse vertebra study. Some test for detecting shape differences are gotten 
and the models are discriminated by the use of a dimension criterion such as the modified 
BIC criterion. 

2 QR Size-and-shape distribution 

Lemma 2.1. Let Y : (N — 1) X K, then there exists a T : (N — I) x n lower triangular 
matrix with ta > 0, i = 1, . . . , min(n, K — 1), and H 6 V Uj k such that Y = TH and 

n 

(dY) = n^- i (dT)(HdH') (2) 



r os sm/2 / A ) 

/ (tr AH) 2t (HdH') = ±1^ £ ^j-C^AA'), 



(3) 
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where C K (B) are the zonal polynomials of B corresponding to the partition k = (ti, . . . t a ) 
of t, with J2T=i ti = t] and (a) K = E[;=i( a - (J - l)/2) ti , (a) t = a(a + 1) • • • (a + 1 - 1), are 
the generalized hypergeometric coefficients and T s (a) — 7r s ( s_1 ^ 4 Ylj =1 T(a — (j — l)/2) is 
the multivariate Gamma fu nction. 

Proof. It is follows from Jamesl (1964), eq. (22) and Muirheadl (1982), lemma 9.5.3, p. 
397. □ 



Theorem 2.1. The QR reflection size- and- shape is 

n 

2 n n nK/2TT t K-i 

_L 1. 11 oo 



/t(T) 



r„ [\k] \n K '* t 



fe (2t) [tr(S~ 1 TT / + »)] ^(flS-^TQ 



t = K 



AT 



where ft = E _1 /i.0 _1 /x', C K (B) are the zonal polynomials of B corresponding to the 
partition k = (ti,... ,t a ) of i, with Ym=i U = t arL d h^\v) is the j-th derivative of h with 
respect to v. 

Proof. The density of Y is 

/y(Y) ! 



|£|*7 2 
1 

|E|*/ 2 

1 

ISF/2 



tr XT^Y - m©- x / 2 )(Y - /I©" 1 / 2 )' 
tr (iT X YY' + S V©" V - 2S" 1 Y©- 1/ V 
tr(S" 1 YY' + n) - 2trST 1 Y0- 1/ V 



with VI = S 1 /x© Taking Y = TH and using Lemma [2~T1 the joint density of H and 

T is 



/h,t(H,T) 



ik- 



tr (E^TT' + ft) - 2tr0- 1/ VXT 1 TH 



Assuming that h(-) can be expanded in power series, see iFang and Zhang I <| 1990T ) . i.e. 



h(a + v)=J2- 



t=o 



Thus 



/h,t(H, T) = [ tr (S- J TT' + n)] [tr (-20" 1 / V'S^TH 



Now, for integration on H G V n ^K, we note that it is zero when t is odd (see Theorem II, 
p. 876 in James (1961) or eqs. (44)-(46) in James (1964)). The marginal density of T is 
expressed as 



/t(T) 



oo 1 



n K/2 U (20! 
J [tr (-20- 1/ VE~ 1 Th) 2 '(HdH'). 
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So, by Lemma [22] and recalling that C\(aA) = a t C\(A), for a constant, we have 



tr ( -20- 1/2 u'S _1 TH 



(HdH') 



r > n 7r nK/2 (I) , 

rl = Y y^-C K f40- 1 /2 /x 'S- 1 TT'S-V©" 1/2 ) 

^ [\ K ] k (l K ) K ^ } 

nn_nif/2 (I) At 

nn nK/2 (1) At 

2 * n V )f&L<7„ (ns- 1 TT / ) . 



From lMuirheadl (|l982h . p.21 ^7t|^ = {k n ) k then Tff = and ^f^T = H> in our 
case k = t, and the result follows. □ 



rf*+a) _ 



Alternatively, the si ze-and-shape density (1H) can be ob t ained as a particular case of the 
singular case studied in iDi'az-Garci'a and Gonzalez-Farias I (|2005l ) . 



3 QR Shape distribution 

Now, observe that for T : (N — 1) x n, n = min(N — 1,K), the matrix T contains (N — 
1)K — nK + n(n + l)/2 non null QR rectangular coordinates (iy ^ 0). Let vechT a vector 
consisting of the no null elements of T, taken column by column. Then the QR shape matrix 
W can be written as 

1 



vechW = - vcchT, r = ||T|| = VtrT'T = ||Y| 



then by Theorem 2.1.3, p.55 of iMuirheadl (|l982t ), 

m / rn \ 

(d vech T) = r m J[ sin m ~ l 9 t I /\ ddi A dr, 



»=i 



with m = (N - l)K - nK + n(n + l)/2 - 1. Denoting u = (0 1; . . .,6m)' and J(u) 

(dT) = r m J(u) dfij A dr. 
Theorem 3.1. The QR reflection shape density is 

n 

2« 7r "^JJ^-V(u) ^ 



/w(W) 



EE 



c re (ns- 1 ww') 



x / r M+2 *- 1 / l ( 2t )[r 2 trS- 1 WW'+trn](dr 
Jo 



(5) 



where M = (N — l)K. 
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Proof. The density of T is 

n 

T\t Ki 

2 n n nK/2 11 » ~ /t( 2t ) [tr(E- lr TT' + SI)] C^ffE^TT') 



r„[x/2] |E|*/ 2 ^ t\ {\k) k ■ 

Making the change of variables W(u) = T/r, the joint density function of r and u is 

n 

2"tt"^/ 2 f =1 f \- [tr^E^WW + SI)] 

/,.w(r,W) - Tn[K/2] |s|K/2 t , 

c,(r»ns-W) 

(PL ( } - 

Now, note that 

• c K (r ! ns- 1 ww') = r 2 *c re (ns- 1 ww / ). 



J. — 1 £= 1 2—1 4 — 1 

• /i( 2 *)[tr(r 2 £- x WW' + n)] =/i( 2 *)[r 2 trS- 1 WW'+trf2]. 

Collecting powers of r by r nK-^±^+2t+ m = r M+2t-i ^ M = (at _ \)K, the marginal of W 
is: 

n 

/w(w) - r.. rlTFl i vi at/2 — 



r n [i^]|s|^ ^ i!(ix) K 

poo 

/ r M+2t - 1 /i( 2t )[r 2 trS- 1 WW'+trf2](rfr). □ 
Jo 



When £ = cr 2 I, then SI = fiQ^n'/a 2 , |S| K / 2 = cr M , and r^rS^WW' = r 2 /cr 2 , 
because tr WW = 1, thus Theorem 13. 1 1 becomes 

Corollary 3.1. The QR reflection shape density is 

n 

/w(w) ~ r n m a* hr^mr~ 

x / r M+2t -^( 2t )[r 2 /a 2 +trn](dr); (6) 



a nd, if the gaussian model is considered with = 1, the resulting density corresponds 
with Goodall an d Mardia I (119931 Theorem 2), see section [5] below. 



4 Distributions excluding reflections 



From subsection 2.1 of Goodall and Ma rdia ] (Il993lh we can derive the QR size-and-shape 



and QR shape densities excluding reflection: 



G 



• If n < K, then Theorems 12.11 13.11 stand for the corresponding T = and W — 
W NR excluding reflection densities. 

• When N - 1 > K and p < K, the QR size-and-shape density for T = is gj) 
divided by 2, where tu > 0, for i — 1,...,K — 1 and txx is unrestricted. When 
N — 1 < K (0| stands, since txx is not present. 

• When N - 1 > K and p < K, the QR shape density for W = W NR is © divided by 
2, and wkk is unrestricted. When N — 1 < K §5§ holds, since wkk is not present. 

• The preceding results also hold when rank/x = K and rankT < K, and event with 
probability zero. 

• However, if p — K, th e excluding reflection d ensiti es do not follow the above rule. F or 
the gaussian case, see iGoodall and Mardia I (|l993l) and lGoodall and Mardial (|l99lh . 

5 Central Case 

The central case of the preceding sections can be derived easily. 

Corollary 5.1. The central QR reflection size-and-shape density is given by 

2"7T^ 

/t(T) = — ^[trE-iTT']. 

r„[iif]|E|f 

Proof. It is straightforward from Theorem just take /x = and recall that M°)[tr •] = 
h[tv}. □ 

And: 

Corollary 5.2. The central QR reflection shape density is given by 

n 

/w(W)= nL^, / rM-'hlrHvV-'WW'Kdr). 

Proof. Just take /x = and &W)[tr ■] = h[tv •] in Theorem [37T1 □ 

The corresponding central excluding reflection densities follows according to Section [U 
Observe that it is possible to obtain an invariant central shape density, i.e. the density 

does not depend on function h(-) Let h be the density generator of Y ~ £n-i,k(0)I® I> h), 

i.e. 

/ Y (Y)=/i(trYY'), 
then bv lFang and Zhang I (|l990t ). p.102, eq. 3.2.6, 



\t»-m- V)dr= niN-i)K/2 
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So, if 8= (trS- 1 WW') 1/2 r, ds = (tr5]- 1 WW') 1/2 (dr), then 

/ r Jlf - 1 ft[r 2 trS- 1 WW'](dr) 
Jo 



M-l , 
fc(s )- 



(trS-iWW) 1 / 2 / '(trS^WW) 1 / 2 

/■OC 

= (trE^WW') - ^ 2 / s M - 1 h{s 2 )ds 

Thus: 

Corollary 5.3. When fi = i/ie Qi? reflection shape density is invariant under the elliptical 
family and it is given by 

1 „ [ 5 AJ 2 . =1 

As in the noncentral case, if £ = ct 2 I, then |S|t = <r M and (trE^WW')"* = er M , 
thus: 

Corollary 5.4. Mien fi = and S = er 2 I £/ie Qi? reflection shape density is invariant 
under the elliptical family and it is given by 

/» = r prW J ( U )I1 W ^ ■ 

1 " L 2 ^ J j=l 

6 Some particular models 

Finally, we give explicit shapes densities for some elliptical models. 
The Kotz type I model is given by 



RT -1 + ^ 

Ky) = - 
r 



So, the corresponding fc-th derivative is 



dy k 



m—l 



JJ(r-l-i) 



i=0 



(-Ry)- 



It is of interest the gaussian case, i.e. when r = 1 and R = |, here the derivation is 
straightforward from the general density. 

K(N-l) 

The required derivative follows easily, it is, h( k \y) = r k{ n_ 1) (— R) k cxp(— Ry) and 



/>oo 

/ r M+2t - 1 h^[r 2 tri:- 1 WW' + tvn}dr 

Jo 



J?* _ M_ t / M \ 

— ^exp(-i?trO)(trS- 1 WW) 2 T — +t . 
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2 n n nK / 2 J(u) Yl 



w K 



Rt _M__ t f M 

x _ W exp(-Rtrn) (trE^WW) 2 T — +i 

27T 2 \ 2 

n 

M - - 

exp(-iZtrn) (trS^WW) 2 J(u)J|w^" 1 



x 



~ r(f +t) ^ c^fls-'ww 1 ) 



So, we have proved that 

Corollary 6.1. The Gaussian QR reflection shape density is 



— M. , r 

etr{-±tt} (trS-'WW) 2 J(u) \\ 



/w(W) = 



i=l 



X 



w/iere M = (JV - 1)K. 

The isotropic case of this density was derived by iGoodall and Mardia (|l993l ) , and it is 

M 

obtained from © noting that C K (aB) = a*C«(B) and, if £ = cr 2 I, then (tr S 4 WW') 2 = 
a M , |S|^/ 2 =a M and that 

c^ins^ww) ° K 1 2a 2 1 - ! 



(trS^WW')' / 1 N ! 



Ck \2a 2 

Finally, we propose the result for the Kotz type I model 



Corollary 6.2. The Kotz type I QR reflection shape density is 



Jw(W) = 



EE 



^(ffiS^WW') 



R r-1 T (M) ^ r S -1 WW ') - ( tr H) r " 1 

^ T~i M-nK . Z~Z~ 

r(f + t + u ) n ^( T _i_ s ) 



J^ u!i2«(trn)«r[r-l + f] 



+ E 

m— 1 

oo 

xE 



m— 1 

JJ(r-l-i) 

u-1 



{-R)- m (trfi)~ 



r[r-l-m+f] 
r[f +* + «] m=l(r-l-m-s)\ 



u=0 



u\R u (trny 



where M = (N - l)K. 

Proof. So, the corresponding fc-th derivative follows from 



n(r-i-i) 



i=0 



(-%)" 



and the associating QR reflection shape density can be obtained after some simplification 
as 



2 n TT nK / 2 JJtu«~V(u) 



/w(W) = 



EE 



_ c^ns-'ww') 

X / r M+2 '-W 2 *V trS ~ lwW ' + trn ]W 
Jo 

n 

,ii ^ v c B (flns- 1 ww / ) 

M _ + 

R T-1 T (M) ( tj . S -1 WW ') 2 ( tr n >) 

2 -«+i 7r M ^ ctr(i?fi) 

r[f + nr^^-i-s) 



T-l 



^ «!i2«(trn)«r[r-l + f] 



+ E 

m— 1 



m— 1 



L i=0 



(-R)- m (tvny m 

r[r-l-m+f] 



E 

M=0 



u\R u (tv£l) u 



. □ 
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6.1 Example: Mouse Vertebra 

This classical application is studied in the gaussian case by iDrvden and Mardia" (Il998h . 
Here we consider again the same model and contrasted it, via the modified BIC* criterion ( 
Yang and Yand (120071) ). with two non gaussian models. 

Here we study three models, the gaussian shape (G), and the Kotz (K) model for t = 
2 and t = 3. The Gaussian shape density was given in ([7]) and the remaining shape 
distributions follows by taking r = 2 and r = 3 in theorem 16.21 However, we need to 
simplify binomial series involved in the terms in braces, this can be done straightforwardly 
but tedious by mathematical induction, the results are summarized as follows. 

Namely, the shape density associated to the Kotz model indexed by r = 2, R = | (and 
s = 1) is given by: 



/w(W) 



M 

Y+ t+l 



7T 2 r„ ( 



x < I tr ^ - 2t ] r 



ctr I 

2 J V 7 t=0 



E 



n 



M 

Y + t 



E 



where M = (N — l)K. 

And the corresponding density for the Kotz model r = 3, is obtained as: 



n 



K- 



l J(u) 



/w(W) = 



£ Mf) 
fj/jj, 

2a 2 



■ etr 



2(J 2 J ^ M (M + 2) 



tr 



2t r 



M 



M 



+ t + l 



M 



+ t + 2 



E— 



*'(f) 



where M = (N - l)K. 

In order to decide which the elliptical model is the best one, different criteria have been 
employed fo r the model selection. W e shall consider a modification of the BIC* statis tic as 
discussed in lYang and Yane ( 2007 ). and which was first achieved by Rissanen (1978) in a 
coding theory framework. The modified BIC* is given by: 

BIC* = -2£(ji, a 2 ,h) + n p (log(n + 2) - log 24), 

where £(/i, cr 2 , h) is the maximum of the log- likelihood function, n is the sample size and n p 
is the number of pa rameters to be estimated for ea ch particular sh ape density. 

As proposed by iKass and Rafterv" ( 1995t ) and Rafterv ( 1995 ) , the following selection 
criteria have been employed for the model selection. 

The maximum likelihood estimators for location and scale parameters associated with 
the small and large groups are summarized in the following table: 
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Table 1: Grades of evidence corresponding to values of the BIC* difference. 



BIC* difference 


Evidence 


0-2 


Weak 


2-6 


Positive 


6-10 


Strong 


> 10 


Very strong 



Table 2: The maximum likelihood estimators 



Group 


BIC* 


/ill 


Ml2 


M21 


fJ-22 


H31 


M32 




a 

K : t = 2 
















if :t = 3 
















-403.824 


1.2398 


39.2181 


13.3663 


3.4263 


22.1414 


-1.4618 


Small 


-418.011 


-3.3846 


44.7126 


14.7682 


5.5268 


25.3360 


1.0451 




-307.863 


3.5716 


131.3120 


44.6939 


11.6686 


74.1405 


-4.5674 




199.6375 


16.9915 


-104.1137 


34.6059 


-4.8256 


65.7152 


17.2009 


Large 


206.7321 


-7.2450 


-90.9671 


28.0714 


-11.2058 


58.7553 


0.8674 




151.6613 


-26.5182 


-71.4992 


20.0230 


-15.3962 


47.4424 


-12.6667 













AHl 


A*42 


M51 


A*52 


a 2 


4.0894 


-4.7493 


-27.1075 


-0.7072 


42.8290 


5.2270 


-4.8965 


-30.7140 


-4.1166 


48.6680 


13.7605 


-15.8399 


-90.7380 


-2.7664 


289.4148 


5.1349 


13.5872 


-82.9587 


-12.7549 


346.0959 


7.1216 


10.3519 


-72.3213 


6.4313 


225.3525 


8.0814 


6.7039 


-56.6913 


21.6223 


109.0523 



According to the modified BIC* criterion (see Table [2]), the Kotz model with parameters 
t = 2, R = o and s = 1 is the most appropriate among the three elliptical densities for 
modeling the data. There is a very strong difference b etween the non gauss i an an d the 
classical gaussian model widely detailed and applied by iDrvd en and Mardia | (jl998l) (and 
previous works) in this experiment. 

Let fi 1 and fi 2 De the mea shape of the small and large groups, respectively. We test 
equal mean shape under the best model, and the likelihood ratio (based on —2 log A ss Xxo) 
for the test Hq . = /j, 2 vs H a : ^ /x 2 , provides the p-value 0.3 10 -12 , which means 
that there extremely evidence th at the mean shapes of the tw o groups arc different. This is 
the same conclusion obtained by iDrvden and Mardia I (jl998l) for a gaussian model. 

A final comment, for any elliptical model we can obtain the Q reflection model, however 
a nontrivial problem appears, the 2t-th derivative of the generator model, which can be 
seen as a partition theory problem. For the general case of a Kotz model (s ^ 1), and 
another models like Pear son II and VII, Bessel, Je nsen-logistic, we can use formulae for 
these derivatives given by ICaro-Lopera et al. ( 20091 ). The resulting densities have again a 
form of a generalized series of zonal polynomials which can be c omputed efficiently af t er som e 
modification of existing works for hypergeometric series (see Koev and Edelman ( 2006h ). 
thus the inference over an exact density can be performed, avoiding the use of a ny asymptoti c 
distribution, and the initial transformation avoids the invariant polynomials of Davis I 
which at present seems can not be computable. 
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